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Abstract 

Since the mid-eighties there has been an accumulation of metahic materials whose thermody- 
namic and transport properties differ significantly from those predicted by Fermi liquid theory. 
Examples of these so-called non-Fermi liquids include the strange metal phase of high transition 
temperature cuprates, and heavy fermion systems near a quantum phase transition. We report on a 
class of non- Fermi liquids discovered using gauge/gravity duality. The low energy behavior of these 
non-Fermi liquids is shown to be governed by a nontrivial infrared (IR) fixed point which exhibits 
nonanalytic scaling behavior only in the temporal direction. Within this class we find examples 
whose single-particle spectral function and transport behavior resemble those of strange metals. In 
particular, the contribution from the Fermi surface to the conductivity is inversely proportional to 
the temperature. In our treatment these properties can be understood as being controlled by the 
scaling dimension of the fermion operator in the emergent IR fixed point. 
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I. INTRODUCTION 



During the last ten years, developments in string theory have revealed surprising and 
profound connections between gravity and many-body systems, resulting in the emergence 
of a new paradigm for strongly coupled many-body systems. The anti de-Sitter/Conformal 
Field Theory (AdS/CFT) correspondence l|] relates a gravity theory in a weakly curved 
{d + l)-dimensional anti-de Sitter (AdS^+i) spacetime to a strongly-coupled d- dimensional 
quantum field theory living on its boundary. This maps questions about complicated many- 
body phenomena at strong coupling to solvable single- or few-body classical problems in a 
curved geometry. Black holes in this geometry now play a surprising and universal role in 
characterizing the strongly coupled dynamics of the boundary theory at finite temperature 
and density, a development anticipated by the discovery of Hawking and Bekenstein in the 
1970s that black holes are intrinsically thermodynamic objects. For example, important dy- 
namical insight into the thermodynamics ^ and transport behavior [sl of strongly correlated 
systems has been obtained from simple geometric aspects of black hole spacetimes. 

Very recently, this apparatus has been brought to bear on the problem of fermions near 
quantum criticality j4-3|. The basic strategy is to perform angle- resolved photoemission 
(ARPES) thought-experiments on a black hole. The ground state of a class of strongly- 
coupled many-body systems is described by such a charged black hole. The fermionic re- 
sponse, which is proportional to the ARPES intensity, may be computed by studying the 
scattering of Dirac particles off of this black hole. See fig. [H for further explanation. In 
particular, by exploring different regions in parameter space, both Fermi liquid-like 6| and 
non- Fermi liquid behavior were discovered, establishing the black hole as an impor- 

tant new tool for addressing outstanding questions related to interacting fermions at finite 
density. 

A prime example of a theoretical challenge to which such a tool may be applied is the 
strange metal phase of the high Tc cuprates. This is a funnel-shaped region in the phase 
diagram emanating from optimal doping at T = 0, and its understanding is believed to 
be essential for deciphering the mechanism for high superconductivity. The anomalous 
behavior of the strange metal — perhaps most prominently the simple and robust linear 
temperature dependence of the resistivity — has resisted a satisfactory theoretical explanation 
for more than 20 years. In particular it implies that the low-energy excitations near the Fermi 
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FIG. 1: An illustration of the holographic approach to understanding non- Fermi liquids. The 
physics of a strongly-coupled many-body theory is exactly dual to that of a gravitational theory in 
a curved spacetime with one extra dimension. A finite density of charge in the many-body problem 
corresponds in the gravity theory to a black hole with an extra emergent conformal symmetry near 
its horizon: these correspond to the symmetries of two-dimensional Anti de-Sitter space. This 
space (when projected onto a two-dimensional plane) forms the basis of the M.C. Escher print 
Circle Limit IV, which is thus used here to represent the horizon of the black hole. Understanding 
fermionic response in this system corresponds to scattering bulk fermions off the black hole and 
shows excitations characteristic of Fermi or non-Fermi liquids. 



surface are not Fermi-liquid-like quasiparticles (8|, |9|, and suggests that the system may be 



governed by a quantum critical point at low energies Similar anomalous behavior 



has also been observed in heavy fermion systems near a quantum phase transition |18N20|. 
where quantum critical behavior is more clearly established. The anomalous behavior of 
a strange metal can be characterized by a phenomenological model called the "marginal 



Fermi liquid" (MFL) js], which assumes that the spin and charge excitation spectra of 
the system are momentum-independent (or have weak momentum dependence) and have a 
specific scale-invariant form.^ 

In this paper we report on a class of non- Fermi liquids discovered using the AdS/CFT 



^ A microscopic model which produces this critical fluctuation spectrum has been proposed recently in [2l| . 
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correspondence. The low energy behavior of these non-Fermi hquids is shown to be governed 
by a nontrivial infrared (IR) fixed point which exhibits nonanalytic scahng behavior only in 
the temporal direction. This fixed point manifests itself in our dual gravity description as 
a region containing a two-dimensional anti-de Sitter spacetime. The scaling behavior that 
arises from it is of the form advocated for the cuprates , and used as an input in 



the study of heavy fermion criticality [22|. Within this class of non- Fermi liquids we find 
examples whose single-particle spectral function and transport behavior resemble those of 
strange metals. In particular, the resistivity is proportional to the temperature. In our 
treatment, these properties arise from a full quantum-mechanical treatment of the system, 
and can be understood as being controlled by the scaling dimension of the fermion operator 
in the emergent IR fixed point. The fact that the sin gle- particle decay rate and the current 
dissipation have the same temperature dependence [sl llj] can be understood holographically 
in terms of the rate at which Dirac particles fall into the black hole (see fig. |6]). 

These results offer hope that new insight into the quantum critical points thought to 
control the behavior of the cuprates and heavy fermion metals may be gained from a dual 
gravity description. 

II. SET-UP 

Many examples of field theories with gravity duals are now known in various spacetime 
dimensions. Well-studied examples include A/" = 4 Super- Yang-Mills theory in li = 4, and 
ABJM theory in ci = 3 23l-l25l|. These theories essentially consist of elementary bosons 
and fermions interacting with non-Abelian gauge fields. The rank of the gauge group is 
mapped to the Newton constant G^v of the bulk gravity such that Gjy oc the classical 
gravity approximation in the bulk thus corresponds to the large A^ limit in the boundary the- 
ory. In addition to these theories, there also exist vastly many asymptotically-AdS vacua of 



string theory 26|, each of which is believed to give rise to an example of the correspondence, 
though an explicit description of the dual field theory is not known for most vacua. 

With a view towards the cuprates, we restrict our discussion to a (2 + l)-dimensional 
boundary theory. For simplicity, we take the field theory to be one with a conformally 
invariant vacuum, which amounts to working with gravity in an asymptotic AdS4 spacetime. 
This conformal symmetry will not play a role in our results below as it will be broken by 
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putting the system at a finite density. 

To set up the gravity description of a system with a finite density of fermions, we consider 
a boundary theory with a f/(l) global symmetry and put the system at a finite chemical 
potential /i for the U{1) charge. This can be described in the bulk by a charged black hole 



whose horizon is topologically in AdS4 spacetime 27|. The conserved current of the 
boundary global U{1) is mapped to a bulk U{1) gauge field Am- The black hole is charged 
under this gauge field, resulting in a nonzero classical background for the electrostatic po- 
tential At{r). For definiteness we take the charge of the black hole to be positive. The zero 
temperature limit of the system is described by taking the zero temperature limit of the 
black hole, which is an extremal charged black hole. In this paper we will be interested in 
zero temperature or more generally T ^ /i. See Appendix |X] for the black hole metric and 
our normalization of gauge field in gravity action. 

Now consider a bulk fermionic field ip. By the general principles of AdS/CFT, this field 
is dual to some fermionic composite operator O in the boundary system. The charge q of 

under the U{1) gauge field can be identified with global f/(l) charge q of O, and the 
conformal dimension A oi O is related to the mass m of ^ by 

A = - + mR (1) 

where R is the AdS curvature radius. When [mRY < y, quanta for if) can be pair produced 



near the horizon 28|. The negative-charged particle in a pair will fall into the black hole, 
while the positive-charged one moves to the boundary of the spacetime. However it cannot 
escape, as the curvature of AdS pulls all matter towards its center, and thus the particle will 
eventually fall back towards the black hole. It then has some probability of falling into the 
black hole or being scattered back toward the boundary. This process will eventually reach 
an equilibrium with a positively charged gas of quanta hovering outside the horizon. See 
figure [21 In this way we set up a finite density of fermions in the bulk: this is dual to a finite 
density of excitations of the fermionic operator O in the boundary theory. 

The charge density carried by the black hole is given by the classical geometry and 
background fields, giving rise to a boundary theory density of order po ~ 0{G'^) ~ 0{N'^). 
In contrast the density for ip quanta is produced from quantum pair production, giving rise 
to a much smaller density of fermions of order pp ~ 0{1) ~ 0{N^). Thus in the large 
N limit, we will be studying a small part of a large system, with the background 0{N'^) 
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FIG. 2: The geometry of an extremal AdS^ charged black hole and the Fermi gas hovering outside 
it. The horizontal line denotes the boundary spacetime, which has d = 2 + 1 dimensions, and the 
vertical axis denotes the radial direction r of the black hole, which is the direction extra to the 
boundary spacetime. The boundary lies at r = oo. The horizon here is a plane, M?. As will be 
discussed in section IIII Al the coordinate r can be interpreted as the energy scale of the boundary 
theory. The black hole spacetime asymptotes to that of AdS4 near the boundary and factorizes 
into AdS2 x near the horizon, with AdS2 including the r, t directions, and M? comprised by the 
spatial X, y directions. 

charge density essentially providing a bath for the much smaller 0(1) fermionic system we 
are interested in. Properties of this finite density fermionic system at strong coupling can be 
found from gravity.^ Here we report results on single-particle spectral function and charge 
transport: 

1. Fermi surface and low energy excitations 

Information regarding the fermionic system can be obtained from the single-particle 
fermion retarded propagator GR{u,k) = (0{u,k)0{—u,—k)\ for O, and the 

\ / retarded 

associated spectral function A{uj, k) = ^lmGR{uj, k). In a physical realization of this 
system A{u,k) could be measured experimentally by Angular Resolved Photoemission 
Spectroscopy (ARPES). In particular the presence of a Fermi surface and its low energy 
excitations manifests itself in the appearance of quasiparticle-like peaks near some shell 
in momentum space. 

The single-particle retarded function Gr{u, k) for O in this strongly coupled system 
can be obtained from the propagator for the bulk field ip with endpoints at the bound- 



For other work on this system, see [29|. 



ary 30|, ISlJ, as shown in figure [3l This amounts to solving the Dirac equation for i{j 
in the charged AdS black hole geometry, which takes the form 



(2) 



where T""^ are gamma matrices in the tangent frame, UabM is the spin connection for 
the metric, and there is a nontrivial gauge field Am due to the charge of the black 
hole. In the full gravity theory, there are also interaction terms for in (jj]), but they 
give rise to higher order corrections in l/N"^ and can be neglected in the large limit. 
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FIG. 3: The two-point function Gii{xi,X2) for a fermionic operator O in the boundary theory 
can be obtained from the propagator for the bulk field ip with endpoints at the boundary. The 
dissipative part of the correlation function, i.e. the spectral function A{(jj,k) = ImG/j(a;, /c), is 
controlled by the rate at which a quantum falls into the black hole, as indicated by dotted lines 
in the figure. In particular, in the low frequency limit, this rate is controlled by the geometry of 
the near horizon AdS2 region. 



2. Charge transport 

The optical conductivity of the fermionic system can be obtained from the Kubo 
formula 

""^"^^ " ^('^-(^)'^-(-^))rctarded (3) 

where is the current density for the global U{1) in x direction at zero spatial 
momentum. The right hand side of ([3]) can be computed on the gravity side from 
the propagator of the gauge field A^ with endpoints on the boundary, as in Fig. HI 
In a 1/A^^ expansion of the conductivity, the leading contribution — of 0{N^) — comes 
from the background black hole geometry. This reflects the presence of the charged 
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bath, and not the fermionic system in which we are interested. As discussed earher, 
these fermions have a density of 0{N^), and will give a contribution to a of order 
0{N^). Thus to isolate their contribution we must perform a one- loop calculation on 
the gravity side as indicated in Fig. |H Higher loop diagrams can be ignored since they 
are suppressed by higher powers in l/N"^. 




+ + + + + + + + 

FIG. 4: Conductivity from gravity. The current-current correlator in can be obtained from the 
propagator of the gauge field Ax with endpoints on the boundary. Wavy lines correspond to gauge 
field propagators and the dark line denotes the bulk propagator for the tp field. The left diagram is 
the tree- level propagator for A^, while the right diagram includes the contribution from a loop of 
quanta. The contribution from the Fermi surface associated with boundary fermionic operator 
O can be extracted from the diagram on the right. 

To conclude this section we note that our gravity analysis of both single-particle spectral 
function and conductivity only involve minimal quadratic couplings of to gravity and the 
f/(l) gauge field, which are governed solely by general covariance and gauge symmetry. That 
is, they concern with a universal sector of essentially all low energy gravity theories. Thus 
our discussion applies to a large class of field theories with an AdS gravity dual with a U{1) 
global symmetry and fermions-the results will depend only on {q,m), and not on the the 
specific nature of the operator nor the details of the theory. In particular, we will treat q 
and m as input parameters, even though in any given theory, only specific values of them 
will arise. The inclusion of other fields in the bulk could change the nature of the ground 
state; we comment further on this issue in the Discussion section. 
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III. RESULTS: NON-FERMI LIQUIDS AND MARGINAL FERMI LIQUID 



A. An Infrared fixed point 

Before describing the results for the spectral function and conductivities, we first discuss a 
feature of the black hole geometry which reveals a crucial dynamical aspect of the boundary 
system. From general principles of AdS/CFT, the coordinate of the extra dimension in the 
bulk, i.e., the radial direction r of the black hole geometry indicated in figure [2|, can be 
associated with the energy scale in the boundary field theory. This connection, along with 
the mapping of isometries of the bulk geometry to boundary theory symmetries, provides 
an organizing principle which associates geometric aspects of the gravity side to quantum 
dynamics of field theory. 

The evolution of geometry from the boundary (at r = oo) toward the interior can be 
interpreted as the renormalization group flow of the boundary theory from high energies 
(UV) to low energies (IR). As r — )■ oo, the black hole metric approaches that of AdS4, which 
has (2 + l)-dimensional conformal symmetries. This reflects that at large frequencies u ^ fi, 
the effects of the finite density become negligible and one recovers the conformal invariance 
of the vacuum. For general r, the presence of the black hole breaks various symmetries of 
AdS4, a reflection of the fact that in the boundary theory at scales w ~ /i, the scaling and 
Lorentz symmetries are broken by finite density. As one moves close to the horizon, one 
finds that the geometry is described by AdS2 x (see caption of figure [2]). The AdS2 factor, 
which involves the time and radial direction of the black hole, has an SL{2, K) isometry, i.e. 
the symmetry group of conformal quantum mechanics, where the scaling symmetry acts in 
the time direction. 

The presence of the near-horizon AdS2 region and associated symmetries indicates that 
at low frequencies {i.e. u <^ fi) the boundary system should develop an enhanced symmetry 
group including scale invariance, and in particular should be controlled by a conformal 
theory, which we will denote the IR CFToi the boundary theory. This IR CFT only involves 
the time direction, with the spatial directions becoming spectators, an aspect which will be 
very important below. We would like to emphasize that the conformal symmetry of this IR 
CFT is not related to the microscopic conformal invariance of the vacuum theory (the UV 
theory) , which at low frequencies is completely broken by the finite charge density. This new 
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infrared conformal symmetry apparently emerges as a consequence of the collective behavior 
of a large number of degrees of freedom. 

It is currently not known how to describe this IR CFT using the conventional language 
of a Lagrangian. Its properties, however, are fully computable from the AdS2 gravity dual. 
For example, the form of the bulk spinor wave equation implies the following. The fermion 
operator with momentum k, 0{k), matches onto an operator Ol^ in the IR CFT with an 
anomalous scaling dimension 6k given by 

Furthermore one can compute the two-point correlation functions of Oj.^ in the IR CFT, 
which we denote by Q^i^)- The emergent conformal invar iance fixes these correlation func- 
tions to be of the form 

Qkiu) = c{k)u^''\ (5) 

where the prefactor c(/c) is complex and analytic in momentum k. It is not fixed by conformal 
invariance and is given explicitly in Appendix |Al Note that due to spherical symmetry both 
Vk and Qk depend only on the modulus k of k. It is important to note that these correlation 
functions Qk{'^) are not the correlation functions of the full theory, which will be described 
in the next section: however, we will see that they directly control the small-frequency 
dependence of the full spectral densities. 

Since we are working with a theory which is scale invariant in the vacuum, all dimensional 
quantities can be expressed in units of the chemical potential /i. In the following, to simplify 
notations we often do not write factors of ii explicitly, which can be reinstated from simple 
dimensional analysis. 



B. One-particle spectral function 



— * 

The retarded function Gji{uj, k) for O in the full theory is found by solving Dirac equa- 



tion ([2]) in the full bulk black hole geometry.^ 
rather involved and are reported elsewhere 



ere we summarize the result. The details are 



Due to spherical symmetry, Gr only depends on = up to a change of basis. We 
first describe its properties at zero temperature. In the low frequency limit u <^ fi, where 
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fi is the chemical potential, it is possible to express Gii{uj,k) in terms of the correlation 
function Qk{i^) for 0{k) in the near horizon AdS2 region plus quantities which depend on 
the full geometry, to which we will refer as UV quantities to distinguish them from objects 
such as (jl]) and that depend only on the structure of the IR CFT. At a generic value of 
k, we find a low frequency expansion 

Gniio, k) = Fo{k) + Fi{k)io + F2{k)gk{io) + ... (6) 

where UV quantities -Fo,i,2(^) are real and analytic functions of k and Qk{^) is given by 
The imaginary part of the correlator and thus the spectral function is controlled by Qk{^)- 
This has a simple geometric interpretation as indicated in figure [31 

For m?R^ < y, the Dirac equation ([2]) can have a static normalizable solution at some 
discrete shell in momentum space, which we call \k\ = kp- Near such a special value of k, 
which can be determined numerically, Gr takes a very different form and its small-frequency 
expansion can be written as 

GR{k,uj) = ^ —— , J:{u,k) = hGkA^) = hc{kF)u"'^^ (7) 

k — kp — —uj — k) 

where Qkpi!^) is given by (ISl) evaluated at A; = kp- The quantities vp.hi and h m. ^ are 
positive constants which are known numerically. The associated spectral function A[uj, k) = 
^lmGR{uj, k) is given by 

A{uj,k) = -- 1^ ^ V ^2 I v2 ' T.{uj,k) = T.i{uj,k)+iT.2{uj,k). (8) 

IT [K — Kp — —UJ — l^i) + 

The non-analytic frequency dependence of the spectral density A{u, k) is again controlled 
by the IR CFT correlation function. 

Equation ([7]) has a pole in the lower frequency plane which approaches oj = Q a.s k kp. 
As a result, the associated spectral function has a quasiparticle-like peak whose location as 
a function of u approaches a; = for — > /c^. The height and width of the peak approach 
infinity and zero respectively. We take this as evidence that the system has a Fermi surface 
at A; = kp. The quantity S(A;,Ci;), which only depends on uj (up to higher order analytic 
corrections in /c — kp), can be interpreted as the self-energy for excitations near the Fermi 
surface, and its imaginary part S2 can be interpreted as the single-particle scattering rate 
near the Fermi surface. See fig. [5] for a density plot of A{uj, k) as a function of k and uj. 
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FIG. 5: Density plot of the spectral function A(u}, k) as a function of uj and k. We have chosen 
units so that ^ = 1 and in all plots m = 0. Left: q = \ with Fermi momentum kp = 0.53 
and Vkp = 0.24 < ^ for which there is no sharp quasiparticle at the Fermi surface and the peak 
dispersion is nonlinear (see equation (jlip ). Middle: q = 1.56 with Fermi momentum kp = 0.952 
and = 0.500 which corresponds to a marginal Fermi liquid. Right: q = 2 with Fermi momentum 
kp = 1.315 and Ukp = 0.73 > ^, for which there are stable quasi-particles at the Fermi surface. 



The two frequency-dependent terms in the downstairs of ([7]) have distinct physical origins. 
The term linear in u comes from an analytic expansion of UV quantities and is real, while 
the self-energy T,{u) comes from the near horizon AdS2 region and is complex, scaling with 
frequency like u'^'^''f . The properties of the quasiparticle-like peak depend on the competition 
between these two terms. The term linear in u in the denominator of can be omitted 
if Ukp < |- For Ukp > I we should still keep the term proportional to uj'^'^^p , since it makes 
the leading contribution to the imaginary part. When Ukp < |, equation ([7]) has a scaling 
form consistent with the scaling hypothesis for a critical Fermi surface discussed recently by 
Senthil 32||, while for v^^ > | it has stable quasi-particles similar to a Fermi liquid. 

More explicitly, when > I? ^he Green's function in equation ([7]) has a pole in the 
lower half complex cu-plane located at 

uj,{k) =uj,{k) -iT{k) (9) 

with 

uj^{k) = vpik - kp) + ■ ■ ■ , ^ (k - kpf'^F-^ ^ 0^ (10) 

uj* yk) 

and the residue of the pole is Z = hivp. The pole represents a quasiparticle which has a 
linear dispersion with vp as the (Fermi) velocity. It becomes stable approaching the Fermi 
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surface, with a non-vanishing spectral weight Z at the Fermi surface itself. In this case the 
standard quasiparticle picture applies even though the scaling exponent of the scattering 
rate is generically different from the w^-dependence of the Landau Fermi liquid. 
For z/fc^ < i, one instead finds that as k is varied toward kp 

uj^{k) {k - kpY , z = - > 1, —y-^ = const (11) 

Also the residue at the pole approaches zero at the Fermi surface 



Z (x{k- kp) '"^^ ^ . (12) 

In this case the imaginary part is always comparable to the real part and the quasiparticle 
is never stable. Also note that the residue of the pole (quasiparticle weight) vanishes at the 
Fermi surface. This is thus an example of a Fermi surface without sharp quasiparticles! 

When z/fc = |, the two frequency- dependent terms in ([7]) become degenerate in the small 
frequency limit. As is often the case for such a degenerate situation, a logarithmic term is 
generated, with the self-energy given by 

'£{u) ^ ciulogu + idiu, ^ = -zi^ (13) 

ci 1 + e ^^ 

where Ci < and di are real constants which can be obtained from (lA4p . Thus in this case 
the single-particle scattering rate is linear in u. While it is still suppressed compared to the 
real part as the Fermi surface is approached, but the suppression is only logarithmic. The 
quasiparticle residue also vanishes logarithmically at the Fermi surface. Remarkably ( fT3|) is 



of the form^ postulated in 



for the "Marginal Fermi Liquid", which fits well with results 



from ARPES experiments on cuprates in the strange metal region [33 1. 

Similar logarithmic terms appear for any z/^^ = ^ ^'^^ example, at z/^^ = 1, 

one finds that the self-energy becomes 

S(a;) ^ cso;^ log u + C2cu^ (14) 



^ Note for marginal Fermi liquids proposed in [8|, |^ is equal to — f for which case there is a particle- hole 
symmetry, i.e. I]2(— w) = E2('^)- In (|13p this happens only for q ~ 0. The particle-hole asymmetry in our 
case for generic q can be attributed to the fact that quantum critical fluctuations are charged. Also note 
in our case (jl3p receives analytic corrections in k in both exponent and prefactors. 
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with C2 real and C2 complex constants. Equation ( !T^ resembles that of a Landau Fermi 
liquid with the imaginary part of the self-energy quadratic in w^, but it has a logarithmic 
term uP' log uj with a real coefficient which is not present for a Landau Fermi liquid. 

At a nonzero temperature T ^ yU, one finds the self-energy T^iiS) in equation ([7]) is 
replaced by 



where Qkp{(^,T) = (27rT)'^'^''F gi is the finite-temperature IR CFT Green's function. 

Qi is a universal scaling function which is given explicitly in Appendix \M It approaches a 
constant for small u/T and gives the correlator at zero temperature in the limit u/T — j- oo. 
At a finite temperature, the pole of (JTj) always has a negative nonzero imaginary part and 
never reaches the real axis. This refiects the familiar fact that Fermi surface gets smeared 
at finite temperature. For the Marginal Fermi Liquid case, z/^^ = i, the self-energy is 
S(ci;) = TiTg2{u) with u = — 4= and (72 (w) is a scaling function given by 



where il) is the digamma function, ci and di are the same constants as in f|T3l) . and ■ ■• 
denotes terms which are real and analytic in uj and T. 

C. DC and optical conductivities 

The contribution of the Fermi surface to the conductivity can be extracted from the one- 
loop diagram in figure HI whose evaluation proceeds similarly to that of a Fermi liquid, except 
for the extra bulk dimension and complications from curved spacetime geometry. There is 
also an important qualitative difference: we are evaluating the diagram in a spacetime 
with an event horizon. For example, the fermion running in the loop can go into the 
horizon as indicated in Fig. [61 In fact, precisely such processes give the leading contribution 
to the dissipative part of the current-current correlator ([3]), and thus the optical and DC 
conductivities. 

Since the system contains a charged background 0{N'^) fiuid, the DC conductivity at 
zero temperature is not really well-defined as there is no mechanism to dissipate away the 
current. In a sense, the system should behave like a perfect conductor. The signature of 
this in the optical conductivity is a delta function at zero frequency with a coefficient of 




(15) 





(16) 
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0{N'^) 351]. We will instead be interested in the leading low temperature contribution from 
the Fermi surface to the conductivity. A well-defined answer for DC conductivity can be 
extracted since at any nonzero temperature of 0{N^), the system contains, along with the 
charged background fluid, a neutral component with density of 0{N'^) which can dissipate 
away the 0{N^) momenta of the fermionic quanta^. 

~ AdS^ 




AdS, X 



FIG. 6: The imaginary part of the current-current correlator ([3|) receives its dominant contribution 
from diagrams in which the fermion loop goes into the horizon. This also gives an intuitive picture 
that the dissipation of current is controlled by the decay of the particles running in the loop, which 
in the bulk occurs by falling into the black hole. 



We find the conductivity (jS]) can be written in terms of boundary theory quantities as 



lUJ 



dk 



2tX 2tX Ui — LO — UJ2 



le 



where f{uj) 



1 



eT+l 



-Aiuji, k) A(a;i, w, k)k{uj2, wi, k) A{uj2, k) 

(17) 

is the Fermi distribution function, Aiu, k) is the single particle spec- 



tral function discussed in last subsection, A is an effective vertex which can be evaluated 
explicitly, and C is a temperature independent overall constant. The explicit derivation and 



evaluation of 



34j are rather complicated and details will appear elsewhere js^ . For an out- 



line of the derivation see Appendix |X1 where an expression for A is also given. Equation ( !T7|1 
can be interpreted in the boundary theory, as seen in Fig. [71 in terms of standard Fermi 
liquid theory where A plays the role of an effective vertex. A is in general a complicated 
function of w, k and T, but in the low temperature limit and near the Fermi surface it be- 
comes a smooth, real function of the modulus k, independent of u and T. In this limit, the 
conductivity is controlled by the single-particle spectral function near the Fermi surface. 



^ The optical conductivity should also contain a piece proportional to 5{ijj) at 0{N^) level, which can be 
interpreted as a correction to the corresponding 0{N'^) piece. The contribution from the Fermi surface 
does not give rise to a delta function at any non-zero temperature. 



15 




A((a,k) 

A{u),k) 
A(-(a,k) 

FIG. 7: Equation (jl7p can be interpreted in terms of a Fermi liquid picture with the exact propa- 
gator for a fermion running in the loop and an effective vertex A. 

We find the DC conductivity, given by a; — ?■ limit of f[T7j) . can be written as 

aDc = a{q,m)T-^'''F (18) 

with a{q, m) a numerical constant depending on the input parameters q and m. In particular, 
for = |, which corresponds to the Marginal Fermi Liquid case, we find the contribution 
to the resistivity from the Fermi surface is linear, as observed for cuprates in the strange 
metal phase and many other non-Fermi liquid materials 



In general, the electrical transport scattering rate can have different temperature and 
frequency dependence from the single-particle scattering rate, because the former emphasizes 
large momentum transfer. Here we find that the scaling exponent in (fT8|) is the same as that 
for S in d?]) and ([H]). On the gravity side, this can be understood intuitively from figure El 
which indicates that the dissipative part of the current-current correlator is controlled by 
the rate of the bulk fermion falling through the horizon, which as described earlier also gives 
the single-particle scattering rate. From the boundary theory point of view, this is possible 
because the single-particle decay rate is controlled by the IR CFT, which is critical only in 
the time direction and does not distinguish between small and large momentum transfer. 

The optical conductivity has the advantage that it can distinguish situations with or 
without quasiparticles, since the existence of a quasiparticle gives rise to a scale which is 
associated with the lifetime of the quasiparticle and is manifest through a transport peak at 
u = 0. We will be interested in the regime u,T fi with u/T arbitrary. More explicitly, 
we find that for u^p < 1/2, where there are no stable quasiparticles, 

a{u) =T-^'''fFi{u/T). (19) 

Here Fi{x) is a universal scaling function approaching a constant as x — )■ 0. For large x it 
falls off as x~'^'^''F which leads to a conductivity of the form, 

a{u) = aiiuy^^^F , T<w</x (20) 
16 



with a a real constant. Note that the phase factor above is fixed by time reversal symmetry 
which requires (J*{oj) = a{—u). The falloff in fj20|) is much slower than the Lorentzian form 
familiar from Drude theory. The behavior (11911201) is consistent with that of a system without 
a scale and with no quasiparticles. 

For > 1/2, there are two regimes. In the first we take u = ujT~'^'^''f fixed in the low 
temperature limit. Then one finds that a can be approximated by a Drude form 



with ul = '^^^, - = 2t;^ImS(u; = 0) oc T^^^F < T (21) 



' ' - — tu ^ 27r ' r 

r 

where C is introduced in (lAllI) in Appendix |A] and vf is the Fermi velocity. As commented 
before the transport scattering rate 1/r is proportional to the single-particle scattering rate 
ImE. When yU ^ a; ^ T, we find the scaling behavior 

a(a;) = ^ + 6(za;)2-..-2 (22) 

OJ 

with b a real constant. The 1/u term gives rise to a term proportional to with a weight 
consistent with that in (12T|) . 

For = the Marginal Fermi liquid, we find 

^iu^) = T-'F,(^^,\og^^ (23) 

where F2 approaches a T-independent constant at small u. Due to time reversal symmetry 
the real part cri^u) of o"(cj) is an even function in u and thus for cj/T < 1, one can again 
approximate cr{uj) by a Drude form with the transport scattering time r oc 7^. For u ^ T 
we find that 

/ N liC I I 1 l + in, 



where ci introduced in (1131) and (1161) is negative. The above expression is consistent with 
that recently obtained for marginal Fermi liquids from summation of the particle-hole ladder 



m 



37|. 



In the strange metal region of the phase diagram for cuprates, it was observed in |l3l. |38| 
that the optical conductivity for Bi-2212 has two different scaling regions: (i) for < T, it 
has a Drude form with ^ ~ T reflecting temperature as the only scale, and (ii) for u > T, 
cr((X') = C{iuj)"'~'^ with 7 ^ 1.33. In region (i) the strange metal appears to be consistent 
with the Marginal Fermi Liquid behavior with Ukp = \- But interestingly in region (ii) the 
strange metal resembles the behavior (122|) for u^p = | > | ^"^^ appears to be consistent 
with that proposed by Anderson j39| based on theories of Luttinger liquids. 
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IV. DISCUSSION 



The physical picture underlying our non- Fermi liquids can be summarized as in figure [8l 
Quantum critical fiuctuations indicated in the figure are described by degrees of freedom 
in the IR CFT and are represented in the gravity side by the near horizon AdS2 region. 
An interesting feature here is that critical fiuctuations remain critical for any momentum, 
similar to that postulated in the Marginal Fermi Liquid description of the cuprates. Here 
it comes from first-principle computations with a full quantum-mechanical treatment of the 
system. 



Quantum 
Critical 



FIG. 8: A cartoon of the underlying physical mechanism of our non- Fermi liquids 



m. 



Non- Fermi liquids also arise from coupling a Fermi liquid to a gapless propagating bosonic 
mode, such as a transverse magnetic excitation of a gauge field, e.g. 4ll, |42|. The forms of 
the fermion Green's functions thus obtained fit into the set of functions we have found in (^^). 
An important difference, however, is that for a gapless boson small momentum scatterings 
are strongly preferred due to its linear dispersion relation. As a result the corresponding 
resistivity grows with temperature as a higher power than the single-particle scattering 
rate |42|. 



The extremal black hole we are studying has a 'residual' zero-temperature entropy. This 
degeneracy is exact in the N ^ oo limit; at finite it could be lifted to a low- lying density 
of states, as in a system with frustration. Nevertheless it is interesting to ask whether the 
IR CFT and its properties are tied to this zero temperature entropy. Such a correlation may 
give a hint about the nature of quantum critical points observed in realistic condensed matter 



systems 



43|. Recently it has been argued in [4^ that when the fermionic backreaction to 



the black hole geometry is included, the black hole horizon is replaced by a geometry with no 
zero-temperature entropy. In the large limit corrections due this replacement are higher 
order in A^^ 401] and qualitative features of our earlier discussion are not affected. 

In a full gravity theory where the extremal charged black hole arises, there could also 
be light charged scalar fields. A mechanism similar to the one described earlier which 
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leads to the bulk Fermi gas in Fig. |2] leads to condensation of those scalar fields, and a 



superconducting state in the boundary theory [45|]. Thus the Fermi surface state and the 
IR CFT we are describing could be hidden behind a superconducting state. In fact, there 
are indications that the existence of such light charged scalars may be generic in string 
theory compactifications 46|. This is not dissimilar to many known condensed matter 
systems where quantum critical points appear to have a tendency to be hidden under some 
superconducting dome. This resemblance may not be an accident and deserves further study. 
The effect of such a condensate on the Fermi surfaces discussed here has been explored 



m 



47 



49|. In such a superconducting state a well defined gapped quasiparticle appears 



for a reasonable coupling between the bulk scalar and spinor 48|. For energies above the 
gap the behavior is as discussed here - an IR CFT is still controlling the physics of the 
quasiparticles despite the presence of the condensate. One can get rid of superconducting 
instability by turning on a magnetic field, which has been discussed recently in 50|-|53|, and 



quantum oscillations in the Fermi surface contribution to the free energy have been observed 



m 



52 



53|. 



To conclude, our models provide a laboratory for studying qualitative features of non- 
Fermi liquids, and offer hope that new insights into the quantum critical points thought to 
control the behavior of the high Tc cuprates and heavy fermion metals may be gained from a 
dual gravity description. For example, the emergent IR fixed point and its consequences for 
single-particle particle and transport may transcend the particular family of gravity models 
that we are studying, and may have applications to real systems. 



Appendix A: Supplementary materials 



1. Gravity action and black hole metric 



The action for a vector field Am coupled to AdS4 gravity can be written as 



2^2 



'7^ + — R^FmnF 



MN 



(Al) 



and R is the curvature radius of AdS. The equations of motion following from (lAlll are 
solved by the geometry of a charged black hole HQ 



(A2) 
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with 



/ = 1 + ^-^, A = Jl-^-], (A3) 



where /i is the chemical potential in dimensionless units. The horizon is at r = 1 and the 
Hawking temperature is T — 



An 



2. Determination of kp and z^^^ dependence on q,m 

A Fermi surface corresponds to a normalizable mode of the Dirac equation in the black 
hole background. Because of the presence of the horizon, this can only occur for static 
solutions, u = 0. Squaring the static Dirac equation and rewriting it as a Schrodinger 
equation one can then reduce the problem of finding a Fermi surface to finding bound 
states in the resulting potential as a function of k. Bound states appear for discrete values 
of = kp. Depending on parameters {q,m) there can be multiple bound states which 
corresponds to the system having multiple Fermi surfaces. For any nonzero frequency the 
Scroodinger potential becomes unbounded at the horizon and the state can tunnel through 
a barrier into the horizon. Estimating the decay rate using WKB one finds the width of 
the state is roughly max(|cjp'^'''-F, T^'^'-'-f) consistent with the self energy of the quasi-particle 
in the field theory. With kp determined, the scaling exponent z/^^ can then be computed 
from dl]). The result is presented in fig. [9l 




-0.4 -0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 



m m 

FIG. 9: Distribution of I'kp in the m — q plane. The two plots correspond to two components for 
a spinor operator in (2 + l)-dimension, which we have been suppressing in the main text. In the 
white region, there is no Fermi surface. There is no quasiparticle in the orange region, v^p < ^- Iii 
the remainder of the parameter space, there is a stable quasiparticle. 
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3. Some long expressions 



The prefactor c{k) in equation ([5]) is given by 



r(-2z/fe) r (i + z/fe - j 



c{k) = e^'-'- ^ ^ '-^^ — , (A4) 



mR 






- l^k 


mR 
V6 




- 7-3- 

''VT2 


+ I'k 



This expression is reminiscent of those for certain two-dimensional CFT. It has been argued 



m 



55 



56[ that the asymptotic symmetry group of the near horizon AdS2 region is generated 



by a single copy of Virasoro algebra with a nontrivial central charge. The scaling function 
gi introduced below equation (fT5|) can be written as 

u k\ r(i + z/fc - 5^ + 



^1 ( ^, - ) = e-^c{k)-^ ^\ (A5) 



where c{k) is given by (IA4I) . 



4. Computation of conductivity 



Here we outline the computation of conductivity f fT7|) . Details will appear in 3J]. Our 
strategy is to first write down an integral expression for the two-point current correlation 
function G^^\iuji) in Euclidean signature and then analytically continue to Lorentzian sig- 
nature inside the integral (for ui > 0) 

G^i^{uj) = Gf{zui = u + ze) (A6) 

where G^^\u) denotes the Lorentzian current correlation function on the right hand side 
of ([3]). The important steps in obtaining G^^\iuJi) are: 

1. In a charged black hole geometry, fluctuations of the U{1) vector field Ar^ which is dual 
to the boundary current Jx mixes with metric fluctuations htx- Thus the wiggly lines 
in fig. m contain the boundary-to-bulk propagators for both vector and metric fields. 

2. In addition to cubic couplings (schematically) ipipA and ipiph indicated in fig. HI the 
couplings of internal fermions to vector and metric fields also include quartic cou- 
plings like h'^ipilj and hAipip, which give rise to seagull-like diagrams. It can be shown 
by explicit calculation that the seguU-like diagrams are always subleading. 
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3. Bulk-to-bulk spinor propagator DE{ri,r2;iuJm,k) for ip as indicated in solid lines in 
fig. m can be written using its spectral decomposition 

DE{ri,r2;tu)m,k) = / , (AT) 

J ZTT lUJm — ^ 

where p{ri, oj, k) is the bulk spectral function for ip, where we have suppressed bulk 
spinor indices. The bulk spectral function p can in turn be related to the boundary 
spectral function Aa/sioJ, k) (jH]) as (we have restored the boundary spinor indices a, 

p(ri, rs; u, k) = V«(n; k) v4„^(a;, k) il)p{r2] u), k) . (A8) 

ij)ct is a normalizable Lorentzian bulk wave function for ip. Again the bulk spinor 
indices are suppressed. The boundary spinor indices a, (3 on i/) label independent bulk 
solutions. 

With the above ingredients and analytic continuation ( 1A6[) . evaluating the one-loop diagram 
in fig. m one finds f|T71) with the effective vertex A given by 



Af}^{ui, 0J2, u,k) = J dr^/^ 'ijjp{r] Wi, k) Q{r; uj; k) il)^{r; 0J2, k) (A9) 
with [Kh and are the boundary to bulk propagators for htx and respectively) 

Q(r; w, k) = -I (-ik^Khir; u)T' + ^drKh{r; w)r*^ - iqKA{r; u)tA . (AlO) 



Note that in flTTI) we have suppressed the boundary spinor indices which are restored here. 

The effective vertex A can be evaluated explicitly. In particular the integration over r 
in flA9l) can be separated into A = A//j + Aj/y, the contributions from the near horizon 
AdS2 region (A/ij) and from the rest of the spacetime (Auv)- One finds that in the low 
temperature limit Auv ~ 0(1) while Am ~ j^^i^fc^+i ^^q^j- ^j^g Fermi surface. 

Since the effective vertex A can be approximated as a constant at the low temperature 
limit, the standard formula to express the optical conductivity in terms of self-energy applies, 
i.e., f|T71) can be written as 

_ / fi-'^-^-fi-') ^ _L_ (All, 

where C is a positive constant and the self-energy S is given in f llSp . 
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